Generation of stable ultra-relativistic attosecond electron bunches via the laser 

Wakefield acceleration mechanism 
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In recent experiments ultra-relativistic femtosecond electron bunches were generated by a Laser 
Wakefield Accelerator (LWFA) in different regimes. Here we predict that even attosecond bunches 
can be generated by an LWFA due to the fast betatron phase mixing within a femtosecond electron 
bunch. The attosecond bunches are stable outside the LWFA and can propagate in vacuum many 
tens of centimeter without significant change in their duration. Our calculations show that evidence 
for the formation of attosecond bunches can be found in the spectrum of coherent betatron radiation 
from LWFA's. 

PACS numbers: 41.75.Ht, 41.75.Jv, 52.38. Kd 
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The generation of attosecond electron bunches would 
be of great interest to provide a new and unique tool for 
modern physical research. The potential of such bunches 
ranges from electron microscopy with attosecond resolu- 
tion, to the generation of attosecond X-ray beams, for 
investigating physical, chemical and biological processes 
on the attosecond timescale. 

The shortest bunches to date are available from Laser 
Wakefield Accelerators (LWFA's) ti]. In recent exper- 
iments exploiting different parameter regimes 0, 0, 0j 
H, @, 0, 0,1^ , ultra-relativistic femtosecond (fs) electron 
bunches were generated by an LWFA. Here we predict 
that even attosecond bunches can be generated by an 
LWFA due to the fast betatron phase mixing withing 
a femtosecond electron bunch. The attosecond bunches 
are stable outside the LWFA, which means that they can 
propagate through vacuum over distances of many tens 
of centimeters, without a significant change of duration. 
We also predict how such bunches can be identified in the 
optical spectrum emitted during the acceleration process. 

Several schemes for the generation of attosecond elec- 
tron bunches have been proposed so far. Examples are an 
inverse free-electron-laser process [13], the interaction of 
high intensity laser pulses with overdense plasma 11| , the 
acceleration of electrons with a short tailored laser pulse 
12|. the slicing of an electron bunch with a laser pulse 
13|, the interaction of an ultra-short laser pulse with a 
nanofilm |,14i] . with a wire or a plasma slice 15|, or the 
interaction of a PW radially-polarized laser pulse with a 
sub- micron droplet of a high-Z material [16| . However, 
all these schemes are of limited attractivity, because they 
either require large accelerator structures [lO| or high in- 
tensities in the order of 10^° — 10^^ W cm~^ (normalized 
amphtude, ao, 10-100) [ll|, [13, E, Q E [ll, such as 
available only from rather exclusive (Petawatt) laser sys- 
tems. A strong disadvantage is that these schemes would 
deliver bunches of rather low energy (a few to a few tens 
of MeV's) which would make it difficult to keep the dura- 
tion of such bunches and apply them due to space charge 



effects. In some schemes the bunches would suffer from 
a limited life time (about 10 fe|l4'l) or a limitation of 
charge to well below a pC [l^jEli 16]. 

Our scheme of attosecond bunch generation, presented 
in this Letter, has the advantage that it is based on the 
nowadays well-known technique of LWFA, and that the 
attosecond bunches are accelerated to GeV energies with 
the use of a Gaussian laser pulse with ao ~ 1 from a com- 
mercially available TW laser system. We show that the 
generated attosecond bunches also live long when they 
have left the laser wakefield. 

The working principle of our scheme is based on a fre- 
quency sweep (chirp) of the betatron oscillation of the 
local bunch radius along an fs electron bunch du ring 
laser wakefield acceleration. Betatron oscillations [17| 
are initiated by a mismatch of the electron bunch radius 
with regard to the focusing field, and the local bunch ra- 
dius will oscillate between some minimum and maximum 
value determined by the bunch emittance and the focus- 
ing gradient [l^] • The oscillation frequency of the radius 
is twice the betatron frequency given by LOp — ^■p\fTJli 
where Wp, / and 7 are, respectively, the plasma frequency, 
the focusing gradient and the relativistic factor [l7|. It 
can be expected that such oscillations occur because the 
bunch becomes naturally mismatched, due to the longi- 
tudinal change in the focusing gradient and due to the 
increase of the relativistic factor during acceleration 17 ] . 
Our investigations show that, despite the extremely short 
bunch duration in an LWFA, the different parts of the 
bunch with different longitudinal positions in the direc- 
tion of propagation undergo betatron oscillations at dif- 
ferent frequencies, mainly due to the longitudinal varia- 
tion of the focusing gradient, but also due to the longitu- 
dinal variation of the relativistic factor. This dynamics 
leads to modulation of the radius of the fs bunch on an 
attosecond scale. This way, sharp peaks of attosecond 
duration are formed in the electron density, during the 
acceleration of electrons to ultra-relativistic energies. 

Our calculations over a wide range of experimental 



parameters, such as typical for channel-guided LWFA's, 
clearly indicate that the described attosecond dynamics 
is a very general and intrinsic feature of LWFA's. More- 
over, the signature of attosecond bunches may even be 
present in the radiation spectra of currently performed 
experiments, as is discussed below. 

In order to demonstrate more details of the described 
bunching in an LWFA, we model the dynamics with an 
fs electron bunch that has been generated via some injec- 
tion mechanism (for example, in the bubble regime [5| or 
using a plasma density gradient injection [7|) and where 
further acceleration occurs in a channel-guided LWFA. 
For the latter, one could tink of a two-stage LWFA, where 
the fs bunch is formed in a gas jet and further accelerated 
in a plasma channel [13] ■ To describe the dynamics with 
typical parameters, we assume that the bunch has ini- 
tially an energy of 51 MeV (7 = 100), an energy spread 
of 3%, a full- width- at-half-maximum (FWHM) duration 
of 7.5 fs, a root-mean-square (rms) width of 1.3 /im in 
both transverse {x and y) directions, and a normalized 
transverse emittance of 1.3 //m. The bunch is injected in 
a 48 mm long capillary discharge plasma channel with a 
radius Vch of 61 /im and an on-axis electron concentration 
of 7 X 10^^ cm~'^ corresponding to a plasma wavelength Xp 
of 40 fxin. The radial unperturbed electron density profile 
np{r) has the typical parabohc profile p^]. The Gaussian 
laser pulse that drives the wakefield in this plasma chan- 
nel has a FWHM duration of 35 fs, is focused to a waist 
of 38 /im, and has a peak intensity of 1.7 x 10^^ Wcm^^ 
(flo = 0.9). The laser pulse is mismatched to the plasma 
channel to suppress the oscillations in its intensity dur- 
ing propagation in the channel due to the self-focusing, 
as was done in [20| . The laser pulse is linearly polarized 
and has a central wavelength of 0.8 /im. The calculations 
were carried out using the fully relativistic particle code 
WAKE 21], which includes also the laser pulse dynamics. 

The electron bunch, which initially has a Gaussian den- 
sity distribution, is injected in the first maximum of the 
accelerating field behind the laser pulse. In Fig. [T] the 
bunch is shown, when it has just left the 48 mm long 
plasma channel, where it has been accelerated to 655 
MeV with an energy spread of 10%. 

As can be seen, the fs electron bunch with initially 
a Gaussian distribution is transformed to a few sub- 
bunches of attosecond duration. The duration of the 
sub-bunches is approximately 600 as. This structure is 
formed due to a transverse change in the distribution of 
electrons, caused by the betatron oscillations of the local 
bunch radius. 

An attosecond structure can be formed in other ways, 
for example, via an external injection scheme like injec- 
tion in front of the laser pulse 2^] or injection at an 
angle [13] • Here we will give an example of the attosec- 
ond structure formed for the scheme where the bunch is 
injected in front of the laser pulse. In this calculation a 
250 fs (FWHM) long electron bunch, which is accelerated 
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FIG. 1: Electron density distribution of the accelerated fs 
electron bunch (in arbitrary units) at the exit of a laser wake- 
field accelerator. 



to a kinetic energy of 2.9 MeV with a standard radio- 
frequency photo-cathode linac, is injected into a plasma 
channel. The bunch has a Gaussian distribution in both 
longitudinal and transversal directions and a 0.7% en- 
ergy spread. We assume that the bunch charge is below 
the beam-loading limit, so that the space charge effects 
can be neglected. The electron bunch is focused into the 
plasma channel to an rms radius of 38 /im. The normal- 
ized emittance of the injected bunch is 0.6 /im in both 
transverse coordinates x and y. The plasma channel and 
the laser pulse parameters are the same as in the first 
example. Immediately after the electrons enter into the 
plasma channel, the high intensity laser pulse follows. In 
the plasma, the electrons will be overtaken by the laser 
pulse, because the velocity of the bunch is smaller than 
that of the pulse. In this process, a considerable part of 
the bunch will be trapped, compressed and accelerated 
in the first accelerating region behind the laser pulse [l^] . 

The typical structure of the trapped bunch that is 
formed and propagates in the wakefield can be seen in 
Fig. [21 Here we see the electron bunch as it comes out 
of the 48 mm long plasma channel. In this distance the 
electrons have been accelerated to an energy of 690 MeV 
with an energy spread of 5.5%. The rms radius of the 
accelerated fs bunch is 1.4 /im and its transverse emit- 
tance is 2.4 /im. Again, a finestructure, caused by the 
betatron phase mixing, is formed inside the bunch. Due 
to the laser pulse dynamics, which especially plays a big 
role during the trapping of the bunch, only one strong 
attosecond density peak shows up. The total FWHM 
bunch duration is 5.4 fs. The FWHM duration of the 
high density sub-bunch is approximately 400 as. 

As was mentioned in the introduction, the stability 
of the attosecond bunches in the vacuum, after leaving 
the LWFA, is crucial. Now we will show that the gener- 
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FIG. 2: Electron density distribution of the accelerated fs 
electron bunch (in arbitrary units) where the initial bunch is 
externally injected in front of the laser pulse. 



FIG. 3: Electron density distribution of the accelerated fs 
electron bunch (in arbitrary units), presented in Fig. [TJ after 
propagating 10 cm in vacuum. 



ated attosecond structure will be preserved in the vacuum 
even after propagation of tens of centimeters. To explain 
this, consider a bunch slice with a fixed longitudinal po- 
sition, z, relative to the bunch center. The radius of the 
slice (the local bunch radius), ct, evolves in the vacuum 
according to [17 1 
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where Z}, = /£„ is the characteristic distance on 
which the radius grows or the "Rayleigh" length and 
Zpr is the propagation distance. The Rayleigh length de- 
pends on the spot size in the focus or the waist cr* , the rel- 
ativistic factor 7 and the normalized emittance e„. The 
waist of the bunch slice is given by crj = 2^/7^/1, where 
h — £^/7^crp+crp^ with dp and (7^ ~ dop/dz the bunch ra- 
dius and the divergence at the exit of the plasma. When 
the bunch leaves the wakefield it can still have focusing 
properties, which means that it will focus at a distance 
behind the channel given by 2:* = (jpu'^jh. The local ra- 
dius in the structured bunch coming from the LWFA can 
be approximated by 17 1 
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where ctq, o\ and C are constants and Ab is the distance 
between neighbor sub-bunches. With the local bunch ra- 
dius and its divergence at the exit of plasma given by 
expression dll), equation (P) can be applied to predict 
the radius of each slice in vacuum. Simple analysis re- 
veals that the Rayleigh length, Zh^ changes monotonically 
throughout the fs bunch having minimum (maximum) for 
minimum (maximum) Up (in these cases a'^ = 0). This 



means that the part of the bunch with smaller radius will 
diverge stronger. From this we can conclude that in the 
vacuum there will be a transition where the structure is 
reversed, which means that the parts with largest radius 
become the parts with smallest radius and vice versa. 
Another way to look at this, is by looking at the trans- 
verse momenta of the electrons, which has a maximum 
for the parts with minimum radius and vice versa. The 
reversing process, confirmed by our simulations, typically 
happens in the first few millimeters behind the plasma 
channel. After this process is finished, the attosecond 
bunch structure stays stable during the propagation over 
many tens of centimeters in vacuum. 

The propagation in vacuum, of the bunch depicted in 
Fig. [ll was calculated using the GPT code [2^ . We found 
that for a typical bunch charge of a few tens of pC, space 
charge effects practically do not play a role and can be 
neglected. We calculated what the bunch profile will be 
after 10 cm of propagation in the vacuum. The result- 
ing bunch density is plotted in Fig. [3l The rms bunch 
radius has grown to 145 /^m, but as discussed above, the 
attosecond structure is maintained, with the difference 
that the maxima and minima of the electron density are 
now reversed. 

The formation of attosecond structure is an intrinsic ef- 
fect of an LWFA and is something that might even occur 
in currently running experiments. However measuring 
the density distribution of an LWFA bunch is something 
that can not be easily done. Here we will show that at- 
tosecond structure in a bunch in an LWFA can be identi- 
fied by looking at the emitted coherent betatron radiation 



24| . We have calculated the radiation for the case pre- 
sented in Fig. [T]the same way as we have done in [23|. In 
Fig. [3]we have depicted the angular and spectral distri- 
bution of the radiated energy normalized to e^/47r^c (here 
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FIG. 4: Coherent betatron radiation emitted by an electron 
bunch in a laser wakefield accelerator. The spectral and an- 
gular distribution of the normalized radiation energy is given 
on a logarithmic scale. 



wakefield accelerated electron bunch can acquire a 
finestructure at the attosecond scale. It turns out that 
this is an intrinsic property of the channel-guided LWFA. 
We showed that structure is formed in a second-stage 
LWFA and in the case of external bunch injection in 
front of the laser pulse. From other simulations we per- 
formed, it became clear that an attosecond structure is 
also formed for external injection at an angle [13] ■ When 
the bunch leaves the wakefield and enters the vacuum, the 
attosecond structure is reversed in a few millimeters of 
propagation. The bunch diverges while further propagat- 
ing in the vacuum, but the attosecond structure remains 
present and is stable during a long propagation distance. 
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e is the elementary charge and c is the speed of hght), 
in the logarithmic scale; 10^ particles were used for this 
simulation. Note that the radiation at frequencies lower 
than LOp will be absorbed in the plasma and cannot be 
observed outside the plasma channel. The radiation is 
confined to very small observation angles in the order of 
l/7m [131, where 7m is the average relativistic factor of 
the electrons at the exit of the accelerator. The strong 
radiation at low frequencies seen in Fig. |4] is the coher- 
ent betatron radiation from the fs electron bunch, which 
we studied in [^J]. This radiation scales approximately 
as exp [— (wCTji/c)^] (here uj is the frequency and Cz the 
rms bunch length) that agrees with the form-factor for- 
malism [25], according to which the energy radiated by a 
bunch is proportional to [1 + {N^ — 1) /], where is 
the number of particles in the bunch and / is the form- 
factor. There is also weaker emission that is peaked at 
higher frequencies, which we found from the simulations 
to be attributed to the fine attosecond structure within 
the fs bunch. To verify this, we have calculated analyti- 
cally the form-factor of a structured fs bunch, modelling 
it as a Gaussian bunch with a local radius that satisfies 
expression Then, we found that the emitted radia- 
tion has a second peak at the frequency 

For the electron bunch depicted in Fig. [U according to 
([3]), there should be a peak in the spectrum of the be- 
tatron radiation at a wavelength of 0.77 ^m, which cor- 
responds to Lu^ « Si2LUp. This agrees very well with the 
value predicted by the numerical calculation (see Fig. [1]). 
In conclusion, our calculations predict that a laser 
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